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The de Broglie-Bohm interpretation of quantum mechanics and quantum field theory is general- 
ized in such a way that it describes trajectories of relativistic fermionic particles and antiparticles 
CO ' and provides a causal description of the processes of their creation and destruction. A general 

method of causal interpretation of quantum systems is developed and applied to a causal inter- 
im ! pretation of fermionic quantum field theory represented by c-number valued wave functionals. 
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^ 1 INTRODUCTION 

+—> • 

The de Broglie-Bohm (dBB) interpretation of quantum mechanics (QM) and quantum field 
theory (QFT) ^ElElEJEl was invented with the intention to provide a causal interpretation of 
the conventional QM and QFT. Yet, the current form of the dBB interpretation of relativistic 
fermions still does not achieve this goal. First, the Dirac spinor is interpreted as a guiding 
wave only for one kind of corpuscles, making no difference between particle and antiparticle 
trajectories (13 El [3]. Second, QFT is not a part of the current form of the dBB interpretation of 
fermions. There exists a causal interpretation of fermionic fields in terms of quantum rotators [SI 
E], but this interpretation works only for free fermionic fields. Besides, the causal interpretation 
of fermionic fields and the interpretation in terms of causal trajectories of corpuscles are viewed 
as two mutually incompatible interpretations. A recent approach [§], based on earlier work 
[TU1 ITT] , partially resolves these problems, but still makes no difference between particle and 
antiparticle trajectories. 

A solution of similar problems for bosons have recently been proposed in Ref. ^2]. In that 
paper, particle trajectories are determined by using the concept of particle current ^3E1E>1 an< ^ 
multiparticle wave functions that result from QFT. A causal interpretation of particle creation 
and destruction is given in terms of causally evolving fields and new hidden variables - particle 
effectivities. Both the particle positions and the fields are considered as beables (i.e., elements 
of reality). The statistical predictions of the theory are compatible with those of the standard 
bosonic QFT. However, despite the compatibility, the predictions of the theory are not equivalent 
to those of the standard bosonic QFT. Instead, the theory has predictions on the statistical 
distribution of particle positions, on which, in general, the standard bosonic QFT does not have 
clear predictions. 

In this paper we generalize the results of Ref. ^2] to fermionic fields. It should be considered 
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mainly as a technical generalization. Consequently, the interpretational issues are not discussed 
in detail as they are already extensively discussed in Ref. ^2]. In Sec. [2] we give a causal 
interpretation of the Dirac equation in terms of particle and antiparticle trajectories. In Sec. 01 
we use the conventional interacting fermionic QFT in the Schrodinger picture [161 1171 118j to 
construct multiparticle wave functions and the equations for the corresponding particle and 
antiparticle trajectories. In Sec. to each fermionic state represented by a Grassmann valued 
wave functional we attribute a corresponding c-number valued wave functional. In Sec. El we 
develop a general method of causal interpretation of an arbitrary quantum system in which a 
c-number valued wave function satisfies a Schrodinger equation. In Sec. El we use the results of 
Sees. 21 and [21 to construct a causal interpretation of fermionic QFT and a causal interpretation 
of particle creation and destruction, in a way similar to that for bosonic fields. The conclusions 
are drawn in Sec. [71 

2 CAUSAL INTERPRETATION OF THE DIRAC 
EQUATION 

The Dirac equation in the Minkowski metric ?7u i / = diag(l, —1, —1, —1) is 

(iffy - m)i>(x) = 0, (1) 
where x = (x ■ X - X . X 3 ) = (t,x). The general solution of Q can be written as 

4,(x) =i; (p \x)+^ A \x), (2) 
where the particle and antiparticle parts can be expanded as 

V (P) (*)=]TW*), ^ A \x) = Y / d* k v k (x), (3) 

k k 

respectively. Here u k (v k ) are positive (negative) frequency 4-spinors that, together, form a 
complete orthonormal set of solutions to (^Q). The label k is an abbreviation for the set (k, s), 
where k is a 3- vector related to a Fourier expansion and s = ±1/2 is the spin label. 

The interpretation of negative frequency spinors as antiparticles (i.e., positrons) above is 
usually considered as a part of second quantization, i.e., QFT. Indeed, QFT is discussed in 
the next section. However, by using the Stiickelberg-Feynman theory, such an interpretation is 
meaningful even without explicit introduction of QFT |19| . 

By introducing the quantities 

n(p)( x ,x') = J2Mx)ul(x'), 

k 

n( A )(x,x ! ) = J2v k (x)vl(x% (4) 
k 

we can extract and ip( A > from ip using 

^ p )(x) = J d 3 x'rt p \x,x')4>(x'), 

x,x')ip(x'), (5) 
where t = t'. We introduce the particle and antiparticle currents defined as 

fp=^ p h^ p \ i A) =^ A h^ {A \ (6) 
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respectively, where tp = V^7o- Since ip( p ^ and tj;^ separately satisfy the Dirac equation (fTJl. the 
currents and are separately conserved: d^jjf^ = d^j\f^ = 0. Therefore, we postulate 
that trajectories of particles and antiparticles are given by 

dx (P) _ j(P)(t, X ( p )) d^ _ j^)(t,x( A )) 
respectively, where j = (J ,j 2 ,j 3 )- 



3 PARTICLE AND ANTIPARTICLE TRAJECTORIES IN 
FERMIONIC QFT 

In fermionic QFT, the coefficients b^ and dt in © become anticommuting operators. The oper- 
ators St (dy create particles (antiparticles), while bk and dk annihilate them. In the Schrodinger 
picture, the field operators V>(x) and ?/>'(x) satisfy the anticommutation relations 

{^(x),^,(x')} = <W5 3 (x-x'), (8) 



while other anticommutators vanish. (Here a is a spinor index.) These anticommutation rela- 
tions can be represented by 



^a(x) 
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<fya(x) 



where rj a (x) and f?„(x) are anticommuting Grassmann numbers: 

{^(x),^')} = K(x),77*,(x')} = {VaWX'W} = 0. 



(9) 



(10) 



Next we introduce a complete orthonormal set of spinors u&(x) and Ufe(x). It is convenient 
to choose them to be equal to the spinors Uk{x) and Vk(x) at i = 0, respectively. An arbitrary 
quantum state may be obtained by acting with the creation operators 



d 3 x$(x)u k (x), 
d 3 x4(x)^(x), 



(11) 



on the vacuum |0) = \^>o) represented by 

*o[»7.*7 + ] =A r exp|y d z x j A'jj t (x)0(x,x>(x')J . (12) 

Here fi(x, x') = ^ (x, x') - Q.^ (x, x'), N is a constant chosen such that (^ol^o) = 1 an d the 
scalar product is 

= Jv 2 rj^*[r],rj^'[r],^}. (13) 

In (|13|) . V 2 rj = DrfDrfi and ^* is dual (not simply the complex conjugate jlbl II Tj ) to ty. The 
vacuum is chosen such that b^o = dk^o = 0. A functional ^[77,77'] can be expanded as 

^h^]=J2 c K^K[r],V t ], (14) 



K 
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where the set {^x} is a complete orthonormal set of Grassmann valued functionals. We choose 
this set such that each \P x is proportional to a functional of the form b kl ■ ■ ■ b\ n d/ k , ■ ■ ■ d/ k , 
which means that each ^> k has a definite number np of particles and a definite number ua of 
antiparticles. Therefore, we can also write (|14j) as 

oo 

*M + ] = *n P ,n A [V,rf]- (15) 

np, n i 4=0 

The tilde on ^ np ,n A denotes that these functionals, in contrast with ^ and ^k, do not have 
unit norm. 

Time-dependent states ^[n,rf,t] can be expanded as 

^[n,rj\t] = J2c K {t)<i> K [r),rf] 
K 

oo 

= E *n P ,n A [V,v\t}. (16) 
np,n A =0 

The time dependence of the c- number coefficients cx{t) is governed by the functional Schrodinger 
equation 

H$,ft]y{ V ,rf,t]=id t iI>[r l ,rf,t}. (17) 
Since the Hamiltonian H is a hermitian operator, the norms 

(*(t)\*(t))=J2\cK(t)\ 2 (18) 

K 

of the states do not depend on time. In particular, if H is the free Hamiltonian (i.e. the 
Hamiltonian that generates the second quantized free Dirac equation ©), then the quantities 
|cr-(£)| do not depend on time. This means that the average number of particles and antiparticles 
does not change with time when the interactions are absent. 

Next we introduce the wave function of particles and ha antiparticles, denoted as 

ip,« A = A 1 -b np d 1 -d nA (xi, • • • , x np , yi, . . . , y„ A , t). (19) 

It has np + ua spinor indices b\, ■ ■ ■ , d nA . For free fields, the (unnormalized) wave function can 
be calculated using the Heisenberg picture as 

i>n P ,n A = (0\^(t,x 1 )--4^(t,y nA )\V), (20) 

where ffi p ) and are extracted from ip with the aid of ©• In the general, interacting case, 
the wave function can be calculated using the Schrodinger picture as 

xCw-Cw*^'''' (21) 

Here the phase (po(t) is defined by an expansion of the form of (|16JI : 

U(t)Vo[v,rf} =r (t)e^>W*oM t ]+ E ( 22 ) 

where ro(t) > and U(t) = U[ip,^,t] is the unitary time-evolution operator that satisfies the 
Schrodinger equation ([170 ■ Note that only ^ np ,n A from the expansion (flUj) contributes to (|2*T|) . 
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The current attributed to the i-th corpuscle (particle or antiparticle) in the wave function 

1pnp,n A is 

MO = I I" Pa t tiCle t S ' i ( 23 ) 

I XI Wi)X ^ or antiparticles, v ' 

where we have used abbreviations 

(24) 



Xb 1 -b np d 1 -d nA = (j^dxdi ■ ■ ■ i r yo)dn A d' nA ^i 1 ... bnp d' 1 -d' n . ' ( 25 ) 



and the superscript T denotes the transpose. Eq. ()23|) is the generalization of © for multiparticle 
wave functions. Therefore, the trajectory of the i-th corpuscle guided by the wave function 
tyn P ,n A is given by the generalization of Q) 

^ = JCQ_. ( 26 ) 

* io(i) 

As an example, consider the wave function describing one positron 

^)*(x)=^(^) = (0|^)t(x)|*). (27) 
From (|23)) we find % = 7oVj X = ^ > so (ESt reduces to 

4 A) = l^ygV = V> (A)t 707^ (A) = ^ (A) 7^ (A) - (28) 

Note that the wave function that corresponds to the presence (not the absence) of a positron is 
ipa^* (not ipa^)- In this sense, we were slightly imprecise in Sec.|2]by referring to as the 
antiparticle part. Nevertheless, the expression @ for the antiparticle current is consistent with 
(|28|) . The wave function Q27|) satisfies the same equation of motion (Dirac equation) as the field 
operator ipi (x). Consequently, the current (|28|) is conserved 



d"i A) = 0, (29) 

and naturally interpreted as the positron current. In particular, the quantity jg = ipa ipa 
is naturally interpreted as the positron density, which can be demonstrated as follows. The 
operator of the positron density is JS] 

n^ A \x)=^ a A \x)^(x). (30) 

Let be a free one positron state 

|*) = $>4|0), (31) 

k 

where ct are some coefficients. The corresponding wave function is given by the complex conju- 
gate of 

^ A Hx)=J24v k (x), (32) 

k 
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while the expectation value of the positron density is 

(V\h( A \x)\*)=iPi A >(x)^ a A Hx). 



(33) 



Note that the positron wave function is not a coefficient of a decomposition of a quantum-field 
state l^) in eigenstates of the positron density operator. In the conventional interpretation of 
QFT, this would imply that the positron density cannot be interpreted as the positron probability 
density. Indeed, the conventional interpretation of QFT does not have definite predictions on 
the probabilities of particle and antiparticle positions ^2]. However, similarly to the bosonic 
case |12| . our interpretation of fermions is an extension of the conventional interpretation. The 
extension consists in postulating the particle and antiparticle trajectories. Eq. (|2fijh together 
with the conservation equation ()29|) . provides the consistency of the probabilistic interpretation, 
i.e., the equivariance. If the positrons in a statistical ensemble are initially distributed with 
probabilities equal to i/^i A ^*^i A \ then (|26|) and (|29|) imply that the positrons will be distributed 
according to 'ipa at all times. 

It is also important to emphasize that, similarly to the bosonic case ^2j, our extension 
is compatible with the conventional interpretation in the sense that any definite prediction of 
the conventional interpretation is also a prediction of our interpretation. With our extension, 
we do not deny any standard QFT effect that can be described purely in terms of quantum 
fields without explicit introduction of the concept of particles. In particular, although our 
extension does not attribute a particle trajectory to the vacuum, the standard measurable effects 
of the vacuum are not denied by our extension. The set of all predictions of the conventional 
interpretation is a subset of the set of all predictions of our interpretation. 



4 BOSONIZATION OF FERMIONIC QFT 

So far, we have not introduced any causal interpretation of the fermionic fields tp or r]. Actually, 
these two quantities are not observables, so there is no need for a causal interpretation of them. 
However, we need a causal interpretation of the processes of creation and destruction of particles 
and antiparticles. For bosonic fields, this can be achieved by introducing a new causally evolving 
hidden variable - the effectivity of a particle ^2] - determined by the causal evolution of bosonic 
fields. Unfortunately, one cannot introduce a similar causal theory for the Grassmann fields rj 
and rf, essentially because the quantity ty*[r],rf ,t]ty{r],rf ,t] is Grassmann valued, so, unlike to 
the bosonic case, one cannot interpret this quantity as a probability density. Therefore, in this 
section we formulate another representation of fermionic states, more similar to that of bosonic 
states. In Sec. using the general method developed in Sec. we use this representation to 
postulate a causal interpretation of fermionic QFT and a causal interpretation of the processes 
of creation and destruction of particles and antiparticles, similar to that of bosonic fields. 

Note first that the notion of the scalar product can be generalized in such a way that it may 
be Grassmann valued |18j . which allows us to write ^[r/, r\\ t] = (f],^^^)} and 

1 = J V 2 r]\r],r]^)(ri,ri Ji \. (34) 

Next, for each fermionic state ^k[v^ V ] we introduce the corresponding bosonic state ^xtty^ V 3 ^! 
with </?(x) and ip' (x) being c-number 4-spinors. The bosonic states are obtained in the same 
way as the corresponding fermionic states, by acting with the bosonic creation operators St and 
d) k (that satisfy the corresponding commutation relations) on the bosonic vacuum ^ol'/ 9 ! <P^]- 
(For details on the bosonic functional methods, see e.g. Refs. jSJ EE EE3 HOI)- Note that the set 
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of bosonic states {^k [Pi <P ]} does not contain states obtained by acting with the same creation 
operator on the vacuum more than once (e.g. (&jj 2 |0)) because such states vanish identically in 
the fermionic case. Therefore, the set of states {^xlv 9 ) is a complete orthonormal basis for 
a subspace of the whole space of all well-behaved functionals <&[<p, ip']. On this subspace, the 
unit operator can be represented as 

l = £l*x> (35) 

K 

We can also introduce the quantity 

k 

= (36) 

K 

so we see that the sets {^k[Vi v ]} an d {^k[<P, P^]} are two representations of the same or- 
thonormal basis {I^A')} for the same Hilbert space of fermionic states. In other words, the 
fermionic state can be represented as *ff[tp, <p',t] = (<p, ip'\^(t)}, which, by using (|34|). (j36|) 

and can be expanded as 

K 

oo 

= E yn P ,n A YP,v\A- (37) 
np,riA=0 

By inserting the unit operator 1 = J T> 2 p\ip, <p')(<p, in (|18|). we see that the time-independent 
norm can be written as 

(¥(t)|tf(*)) = J TP<py*[<p,<p\tySl{<p,<J ,t\. (38) 

Therefore, we interpret the quantity 

p[(p,<p\t] = V*[ip,tp\tyB[<p,ip\t] (39) 

as a positive definite probability density for spinors cp and ip' to have space dependence y(x) 
and ^(x), respectively, at time t. 

The Schrodinger equation (|17[) can also be written in the ^-representation as 

H v V[<p,<p\t] = id t V[<p,<p i ,t]. (40) 

The Hamiltonian Hm is defined by its action on wave functionals *ff[ip, (p* , t]. This action is 
determined by 

H v *[<p,<pi,t] = Jv 2 V Jv 2 p'(p,^\ V ,^) 

xH{^y^'^y^'\ti (4i) 

where H = H\ip,\p^\ is the Hamiltonian in (|17|). 

Having constructed the ^-representation of fermionic QFT, we can introduce a causal inter- 
pretation of it. We first develop a general method of causal interpretation in the next section 
and then apply it to fermionic QFT in Sec. H3 
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5 GENERAL METHOD OF CAUSAL INTERPRETATION 



In this section we develop a general method of causal interpretation of a quantum system de- 
scribed by a c-number valued wave function that satisfies a Schrodinger equation. For simplicity, 
we study the case of n degrees of freedom denoted by a real re-dimensional vector <p, but we 
write all equations in a form that can be easily generalized to the case with an infinite number 
of the degrees of freedom. 

The wave function ip((p, t) satisfies the Schrodinger equation 

Hip = id t ip, (42) 

where H is an arbitrary hermitian Hamiltonian (written in the (^-representation). The quantity 
p = ip*ip is the probability density for the variables (p. The corresponding average velocity is 



dim) 

dt 

where u is defined as 
We introduce the source J defined as 



(P<pp(0,t)ii((?,t), (43) 



U = Rei —^T- (44) 

J=|!+V(pre). (45 ) 

For many physically interesting Hamiltonians, the source J vanishes. However, in general, J 
does not need to vanish. An example for which J ^ is provided by the Hamiltonian H v defined 
by (|41|) . We want to find the quantity v(ip,t) that has the property (|4*H|) 



dim) 



(r<pp(<p,t)v(<p,t), (46) 



dt 

but at the same time satisfies the equivariance property 

^ + V(pv) = 0. (47) 

These two properties allow us to postulate a consistent causal interpretation of quantum me- 
chanics in which (p has definite values at each time t, determined by 

^ = v(m- (48) 

In particular, the equivariance Q47|) provides that the statistical distribution of the variables ip 
is given by p for any time t, provided that it is given by p for some initial time to- When J = 
then v = u, which corresponds to the usual dBB interpretation. The aim of this section is to 
generalize this to the general case J ^ 0. 
Let us write the velocity v in the form 

v = u + p- 1 S, (49) 

where £(<p,t) is the quantity that needs to be determined. From Q47JI . we see that £ must be a 
solution of the equation 

V£ = -J. (50) 



S 



Let E(if,t) be some particular solution of (|5Uj). Then 



£(<p,t) = e{t) + E(<p,t) (51) 

is also a solution for an arbitrary (^-independent function e(t). Comparing (|46|) with l|43|l. we 
see that we must require / d n ip£ = 0. This fixes the function e(t) to be 

e(t) = -V- 1 f d n pE{<p,t), (52) 



where V = J d n ip. It remains to choose the particular solution E. We choose it such that E = 
when J = 0. Eqs. (|52[) and (|51[) then imply that £ = when J = 0, i.e. that v = u when J = 0. 
There is still some arbitrariness in choosing the particular solution E, so we proceed in a way 
that seems to be the simplest one. We take E to be a gradient of a scalar function 

E = V$, (53) 

so $ satisfies the Poisson equation 

V 2 <£ = -J. (54) 

It is convenient to solve (|54|) by using the Green-function method. By taking the solution to 
have the form 

$(M= ! dr<p'Gfa<j?)j{0,t), (55) 



we see from (|54|) that G satisfies 

V 2 G(<p,0') = -5 n (<p-0). (56) 
The solution of ()56|) is most easily expressed as a Fourier transform 

^ f d n k 

G( ^ ) = i(^^ (57) 

To eliminate the factor l/(2-7r) n , we introduce a new integration variable x = k/2ir. In this way, 
(f55|) becomes 

*(<p,t) =Jd n xJ <r<p' J(/,t), (58) 

which is written in the form appropriate for a generalization to an infinite number of the degrees 
of freedom. 

The formalism above can also be easily generalized to the case in which the integration 
measure in (|43[) is not simply d n ip, but a more general expression of the form d n <p fj,(<p,t). In 
some cases, [i is a Jacobian that can be eliminated by a redefinition of the integration variables 
(p. If this is not the case, [i can always be absorbed into a redefinition of the density p — > /ip. 

There is one additional tacit assumption in our general analysis that does not need to be 
satisfied in an even more general situation: the components ipi of the vector (p are assumed 
to be continuous parameters. If these parameters were discrete instead of continuous, then 
Eq. (|43|) would contain a summation over allowed values of <p instead of the integration over 
them. Examples of discrete variables are the z-component of spin and the momentum of a 
particle confined in a finite volume. There exists a consistent causal interpretation for these 
two examples 5,, but it is not clear how to introduce a general systematic method of causal 
interpretation of discrete quantum variables. Formally, one might reduce the summation to an 
integration by introducing the integration measure in which the function p, is essentially a sum of 
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5-functions. However, owing to the singular nature of 5-functions, it is not clear that this would 
lead to a mathematically well-defined causal theory (see e.g. Ref. [2*T]). A recent approach [22] 
suggests a solution of this problem based on smearing of the (^-functions, but this approach does 
not seem to be unique. Therefore, the formulation of a general method of causal interpretation 
of discrete quantum variables is still a challenge. Fortunately, we do not need it in this paper. 



6 CAUSAL INTERPRETATION OF FERMIONIC QFT 

We are now ready to formulate a causal interpretation of fermionic QFT. In Sec.HJwe have been 
working with complex spinors tp and <p* , having complex components cp a and ip*, respectively. 
On the other hand, in Sec. |S] we have developed a general method of causal interpretation for 
real variables. Therefore, in this section we first introduce real spinor components ip\ and (p 2 
defined by <p a = {ip\ + iip 2 )/y/2, so that 

1 <P* a + <Pa 9 • & ~<Pa r, Q s 

*a = -^, Va = t—^- (59) 

Each functional of the complex spinors ip and tp< can also be viewed as a functional of the real 
spinors ip 1 and ip 2 . To simplify the notation, when appropriate, we denote the functionals of the 
form of ^[i^ 1 , ip 2 ,t,x] by a shorter simbol A[x]. 

Following the general method developed in Sec. we introduce the quantity 

where r = 1, 2 and \P = ^[v? 1 , ty? 2 , t]. Next, we introduce the source 

dt ^a(x) V ' 

where p = ^*^/. Introducing the notation 

= W ^«I(x)^(x), (62) 

a,r 

Eq. (|58|) generalizes to 

^\<p 2 ,t] = / V 2 X / V\' e —-- 2 JW\v' 2 A. (63) 

Introducing the quantities 

E r a [x] = ^rr, (64) 

e£(t,x) = -V- 1 J PV^b 1 ,^,^], (65) 

where V = J T> 2 ip, the corresponding velocity is given by 

<[x]=<[x]+p- 1 (e^(t,x)+^[x]). (66) 

Now we introduce hidden variables (/? 1 (t,x) and ip 2 (t,x). The causal evolution of these hidden 
variables is given by 

d<Pg(t,x) 1 2 

-g t = v a [tp ,(p ,t,x], (67) 
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where it is understood that the right-hand side is calculated at <^ r (x) = ip r (t,x). 

In analogy with ordinary bosonic fields ^2], we also introduce the effectivity e np ,n A °f the 
particles guided by the wave function ip np ,n A , given by 

\^ ,o 2 f\ - Ijw^WiiP (en 

Since the physical meaning of the effectivity is discussed in detail in Ref. f° r ordinary 
bosonic fields, here we only note the final results valid also for the bosonic representation of the 
fermionic states above. In an ideal experiment in which the number of particles is measured, 
different ^ np ,n A functionals do not overlap in the (</r, (p 2 ) space. Consequently, the fields (p 1 and 
ip 2 necessarily enter one and only one of the "channels" ^ np ,n A - In this case, the effectivity ()68|) 
of the corresponding particles is equal to 1, while that of all other particles guided by other wave 
functions ip n ' p)n ' A is equal to 0. The effect is the same as if the wave functional VP "collapsed" 
into one of the states \P np nA with a definite number of particles and antiparticles. This provides 
a causal description of the creation of particles guided by ip np ,n A and the destruction of all other 
particles guided by other wave functions ip n ' n ' ■ 

Note that both the corpuscle (particle and antiparticle) positions and the commuting fields 
(p , ip 2 are the preferred hidden variables that describe the elements of reality associated with 
fermions. This would not be possible if the corpuscle positions were described by hermitian op- 
erators that do not commute with tp , ip 2 . However, just as in the relativistic bosonic case |12j . 
there are no hermitian corpuscle-position operators associated with the corpuscle positions. At 
the level of first-quantization formalism, the absence of the hermitian corpuscle-position opera- 
tors is a consequence of the fact that the separate sets {ufc(x)} and {vk(x}} are not complete, so 
one cannot construct a (5-function as a linear combination of the functions of only one set. Sim- 
ilarly, at the level of second-quantization formalism, the are no 1-particle states |^^) with the 
property h^ p \x)\^^) = 5 3 (x)\^f^) (and analogously for 1-antiparticle states). These facts 
are also related to the already mentioned fact that the conventional theory, in general, does not 
give definite predictions on the probabilities of corpuscle positions. In our theory, the corpuscle 
positions correspond to an extension of the conventional theory. As discussed in Sec. EJ the 
consistency of the probabilistic interpretation of particle and antiparticle wave functions rests 
on the fact that the probabilities correspond only to an effective statistical description of the 
underlying deterministic Bohmian trajectories. 

7 CONCLUSION 

In this paper, three new results related to the causal interpretation of fermions have been 
obtained. First, the usual causal interpretation of the Dirac equation is modified such that it 
makes a difference between particle and antiparticle trajectories. Second, fermionic interacting 
QFT is used to construct the multiparticle wave functions that describe particle and antiparticle 
trajectories in systems with many particles. Third, a causal interpretation of fermionic fields 
themselves, based on bosonization of fermionic fields, is constructed and used to formulate a 
causal interpretation of particle creation and destruction. 

Besides this progress in the causal interpretation of fermions, a general method of causal 
interpretation is developed. It can be applied to any quantum system described by a c-number 
valued wave function ip(<p,t) satisfying a Schrodinger equation, provided that the c-number 
parameters ipi are continuous. 
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